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ABSTRACT 



With proper physical mechanisms of energy and momentum input from around the 
centre of a self-gravitating polytropic gas sphere, a central spherical "void" or "cavity" 
or "bubble" of very much less mass contents may emerge and then dynamically expand 
into a variety of surrounding more massive gas envelopes with or without shocks. We 
explore self-similar evolution of a self-gravitating polytropic hydrodynamic flow of 
spherical symmetry with such an expanding "void" embedded around the center. The 
void boundary supporting a massive envelope represents a pressure-balanced contact 
discontinuity where drastic changes in mass density and temperature occur. We obtain 
numerical void solutions that can cross the sonic critical surface either smoothly or 
by shocks. Using the conventional polytropic equation of state, we construct global 
void solutions with shocks travelling into various envelopes including static polytropic 
sphere, outflow, inflow, breeze and contraction types. In the context of supernovae, we 
discuss the possible scenario of separating a central collapsing compact object from an 
outgoing gas envelope with a powerful void in dynamic expansion. Initially, a central 
bubble is carved out by an extremely powerful neutrinosphere. After the escape of 
neutrinos during the decoupling, the strong electromagnetic radiation field and/or 
electron-positron pair plasma continue to drive the cavity expansion. In a self-similar 
dynamic evolution, the pressure across the contact discontinuity decreases with time to 
a negligible level for a sufficiently long lapse and eventually, the gas envelope continues 
to expand by inertia. We describe model cases of polytropic index 7 = 4/3 — e with 
e > and discuss pertinent requirements to justify our proposed scenario. 

Key words: hydrodynamics - ISM: bubbles - ISM: supernova remnants - shock 
waves - stars: winds, outflows - supernovae: general 



1 INTRODUCTION 

Voids of much less density as compared to more massive 
and grossly spherical surroundings are fairly common in as- 
trophysical systems on various spatial and temporal scales. 
They have been found in a wide diversity of settings, in- 
clu ding planetary nebulae (PNe) (e.g. NGC 7662 and NGC 
Lou fc Zhai 2010), bubbles and 



40; Guerrero et al. 



superbubbles (e.g. 



2004 



McCray fc Kafatos|[T987 Korpi et al 



19991 in the interstellar medium (ISM), supernovae (SNe), 
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supernova remnants (SNRs) and so forth. It is conceivable 
that with sustained powerful sources of energy and mo- 
mentum released from around the central region of a self- 
gravitating gas sphere, a considerably rarified void or cav- 
ity or bubble can form and dynamically expand into a more 
massive surrounding envelope. Depending astrophysical con- 
texts, such sustained central sources could be tenuous stel- 
lar winds, magnetized relativistic pulsar winds of mainly 
electron-positron pair plasma, energetic neutrino flux, and 
radiation field of trapped photon gas etc. We advance in 
this paper a theoretical model scenario, which is formulated 
within the framework of self-similar evolution for a conven- 
tional polytropic gas, towards general profiles of spherically 
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symmetric hydrodynamic systems embedded with central 
voids in expansion. As a simplifying approximation, these 
voids are treated as massless quasi-spherical cavities whose 
gravitational fields are negligible when considering the dy- 
namic evolution of the gas shells surrounding the central 
voids. Naturally, there are always some materials inside the 
boundary of any void in reality. They are even indispensable 
for explaining the dynamic evolution of the shell outside in 
aspects other than gravity (e.g., pressure driving, tenuous 
wind driving, extremely relativistic light particles, and pho- 
ton gas etc.). With these qualifications in mind, our analysis 
would indicate that this model is capable and applicable for 
describing a variety of polytropic gaseous astrophysical en- 
velopes in hydrodynamic evolution. Some similar analysis 



along this line can be found in (Lou & Zhai 2009 2010), 
which focus on isothermal cases (i.e. isothermal self-similar 
voids - ISSV) and astrophysical applications to PNe NGC 
40 and NGC 7662 as examples; such central ISSV in PNe 
are powered by hot tenuous stellar winds with shocks. 

Originated from extremely violent processes like super- 
nova explosions, "hot bubbles" or voids filled with intense 
radiation fields are a significant kind of structures in the evo- 
lution of many astrophysical systems including supernova 
remnants (SNRs). One of our main concerns here is on var- 
ious void structures in the self-similar dynamic evolution of 

fl999) noted 



SNe and subsequent SNRs. Matzner & McKee 



the existence of similarities in the evolution of SNe at a 
very early epoch - much earlier than the Sedov stage known 
for its characteristics of self-similar evolution. It would be a 
first approximation to investigate dynamic evolution of SNe, 
whose mass is sufficiently large ( > 8Mq) as compared with 
the mass of the central collapsed compact object (~ 1.5Mq), 
by self-similar hydrodynamic models with central voids. 

Formation of such initial "voids" within one or two hun- 
dred kilometers is plausibly the result of intense neutrino 
flux heating and driving, leading to the emergence of a re- 
bound shock. This Wilson mechanism was put forward and 
elaborated by [Bethe fc Wilson] ( |1985[ ) (and subsequently in 
Bethe 19901, where the authors argued that acceleration 



of mass infall towards the center is inhibited by the neu- 
trino pressure from the so-called neutrinosphere with little 
amount of mass. This mechanism is highlighted in numeri- 



cal simulations of Janka & Hillebrandt ( 1989b I and Janka 



& Miiller (19961 for examples. These simulation works have 
sketched a picture by simulation that stellar materials are 
heated and pushed outwards by the intense neutrino flux 
generated by nuclear processes during the rapid core collapse 
at the centre. As a by-product of rebound shock revitaliza- 
tion, a bubble or cavity will have already been shaped up 
around the center during the epoch before the surrounding 
gas materials become too tenuous (~ 10^ g cm~^ as dis- 
cussed presently in subsection 4.3 1 to trap extremely ener- 



getic neutrinos. Then the envelope shell, consisting of the 
vast majority of stellar materials, is expected to be ejected 
by the revived rebound shock. After a short while, relativis- 
tic neutrinos and outer surrounding massive gas envelope 
becomes completely decoupled. 

We further discuss this physical feasibility regarding two 
factors: energy and momentum transfers from photons to an 
ionized plasma are very effective (compared with that from 
neutrinos to gas); the energy carried by radiative emissions 
of photons is of the same order of magnitude (~ 10'''^ erg) 



as the energy needed to blow a stellar envelope up. Our 
dynamic void solutions show that it could be the photon ra- 
diation field trapped inside the central cavity that drives an 
outer shell to expand after the decoupling and escape of en- 
ergetic neutrinos. As simplifications, we would treat a void 
containing a uniform radiation field instead of a mixture 
of conventional matters and radiation coupled by transport 
equations. The scale of a central void at the very early stage 
( > 100 km, e.g. Bethe|l990 1 makes this uniformity approx- 
imation plausible as the perturbation in radiation field trav- 
els at the speed of light c. Here the concept of radiation field 
is usually generalized as a combination of photons (electro- 
magnetic field) and various products of pair-production such 
as electron-positron pairs, since the temperatures are always 
sufficiently high (at least fceT > 1 MeV) for pair production, 
where fcs is the Boltzmann constant. This radiation-driven 
envelope expansion with self-gravity is then responsible for 
the acceleration and/or deceleration in the expansion of a 
SN; we shall come back to this possibility later in this paper. 

We conduct in this paper a systematic analysis of the 
self-similar hydrodynamic evolution of void surrounded by 
spherical envelopes with various radial structures using the 
similarity transformation. Self-similar hydrodynamics with 
spherical symmetry has been extensively studied as a valu- 
able tool with various approximations and simplifications. 
The research works of Larson! (1969a b|^ Penstonj (1969a'b]), 

Hunter] ([lff7l[l986ll,[Shunm fc Hsu|(|1995|),,Cheva- 

'■- ■'- - - 



lier| (|1997|), |Shu et al.| ([2002[), [Lou"fc Shen[ ( |2004| , |Shen 



& Lou (20041, and Bian & Lou (20051, have investigated 



self-similar hydrodynamics and their astrophysical applica- 
tions with an isothermal equation of state (EoS). Meanwhile, 



those with polytropic EoS have been studied by|Goldreich 


& Weber (19801, Yahfl (1983 


1, Lattimer et al. (19851, 


Suto 


& Silk (19881, Lou & Wang 


20061, Hu & Lou (20081 


and 
3 our 


Lou & Cao ( 2008 1. Some of them are more pertinent t( 



investigation here, e.g. Hunter (19771 obtained a system- 



atic self-similar transformation and obtained the isother- 
mal expansion- wave collapse solution (EWCS), which has 



been substantially generalized; Tsai & Hsu ( 1995 1 intro- 
duced a procedure of treating shocks in an isothermal gas; 



Yahil ( 1983 1 and Suto & Silk ( 1988 1 developed self-similar 



procedure for conventional polytropic gas dynamics; |Lou fc| 
Wang ( 2006 1 and Lou & Cao ( 2008 1 explored the scheme for 



the evolution of conventional and general polytropic gases 
with shocks. Self-similar transformations permit a series of 
central- void solutions ( Lou fc Zhai|2009 2010 for isothermal 
cases and Hu fc Lou|2008 for polytropic voids as examples) 
whose boundaries appear to be very steep in mass density. 
This treatment is again an idealization and the justification 
of doing so is discussed in Appendix |D] We are also inter- 
ested in the construction of self-similar conventional poly- 
tropic shock solutions with expanding central voids. 



We note |Yahil| (|1983| and |Lou fc Cao| (|2008| for their 
analyses of cases whose polytropic indices are 7 — >■ (4/3)"'' 
and 7 = 4/3. For either relativistically hot or degenerate 
materials, it is a very good approximation to describe their 
behaviour by an adiabatic EoS p = Kp^^^ . Here k is a func- 
tion of (r, t) in general polytropic cases satisfying specific 



entropy conservation along streamlines (e.g. Fatuzzo et al 



2004 1. Specification of global constant is a special case, 

i.e., a conventional polytropic gas. An exact 7 = 4/3 case 



is analyzed by Goldreich fc Weber ( 1980 1 for a homologous 
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core collapse solution. In this paper, we emphasize void solu- 
tions with conventional poljrtropic conditions whose 7 takes 
the form of 7 = 4/3 — e with a small e > and discuss the 
physical implication of e. We propose possible situations to 
sustain such void expansions in astrophysical contexts. 

This paper is structured as follows. Section[l]introduces 
background information, the physical problem, and our re- 
search motivation. In Section [2] we show and discuss the 
formulation of self-similar transformation of hydrodynamic 
equations, as well as solution behaviours near sonic criti- 
cal points and shock jump conditions. Section |3] discusses 
the behaviours of solutions near the void boundaries and 
gives various conventional polytropic self-similar void solu- 
tions crossing the sonic critical surface by different manners; 
shock solutions with different kinds of dynamic envelopes are 
presented. We shall explore applications of those models es- 
pecially in the context of SN in Section|4]and some examples 
are presented there. Section [5] contains conclusions and re- 
lated discussions. Appendices A to E offer details of physical 
reasoning and some mathematical derivations. 



2 FORMULATION OF SELF-SIMILAR MODEL 

We first present below the nonlinear Euler hydrodynamic 
partial differential equations (PDEs) with spherical symme- 
try and self-gravity. In spherical polar coordinates in which 
a spatial point is defined by (r, 9, <j}) where r is the radius, 
9 is the polar angle and is the azimuthal angle, we will in- 
troduce self-similar transformation for a conventional poly- 



analysis [e.g. Suto & Silk ( 1988 1] for a conventional poly- 



tropic gas (e.g. Suto fc Silk 1988 Hu fc Lou 2008 1 . For spher- 



ically symmetric hydrodynamics, all the dependent physical 
variables do not vary with angles 9 and (j>. 

2.1 Euler hydrodynamic PDEs, conventional 
polytropic self-similar transformation, 
and asymptotic solutions 

With spherical symmetry and self-gravity, we have the radial 
momentum equation as 

1 dp _ GPI 
p dr ^2 ' 

where we use u to denote the radial flow velocity, t for the 
time, p for the pressure, p for the mass density, M for the 
enclosed mass within radius r at time t, and G — 6.67 x 
10~^ dyn cm^ g~^ for the universal gravitational constant. 
Nonlinear PDE ([T]) is coupled with the mass conservation 



du ^ du 
dt dr 



(1) 



dt Or 



, 



(2) 



or equivalently, in the integral form of 

dM dM „ 
dt dr 

dM 2 
— — = 47rr p . 
dr 

A closed set of nonlinear PDEs is established when we in- 
clude the conventional polytropic EoS, 



(3) 



P = Kp 



(4) 



where 7 is a constant index and k is a global constant. 

A self-similar transformation based on dimensional 



tropic gas is introduced below, namely 

u{r,t) = k^^'^t"^^v{x) , p{r,t) = 



- = W^"^^^^ ' 

(5) 

where k is the sound parameter and n is the scaling expo- 
nent, X is the dimensionless independent similarity variable, 
and depending on x only, v{x), a{x) and m{x) are the di- 
mensionless reduced variables for radial flow velocity u{r, t), 
mass density p(r, t) and enclosed mass M{r, t), respectively. 
Under transformation (|5|, nonlinear PDEs |T]) to Q are 
converted into two coupled nonlinear ordinary differential 
equations (ODEs), viz., 

da _ a 

dx (nx — v)'^ — 7aT^i 



X 

dv 
dx 



{x — v)(nx — v) 



(3n - 2) 



(6) 



v){n - 



l)v + 



{nx - 



(3n - 2) 



27: 



According to Suto & Silk ( 1988 1, the conventional polytropic 
condition for a time-independent k in polytropic relation Q 
simply requires n -4- 7 = 2. 

With n -I- 7 = 2 and the assumption that \v{x)\ and 
a{x) are non-increasing functions of x at very large x, v{x) 
and a{x) have the following asymptotic behaviours^ in the 
regime x^ 1, namely 



a = Ax 



-2/n 



V = Ax 



(n-2)/n 



2(2 - n) 

nA" 
-1) 



(3n-2)_ 



(7) 



(Lou & Shi 2011 in preparation), where A and B referred to 
as mass and speed parameters are two integration constants, 
characterizing asymptotic behaviours at large x. 

The isothermal case of n = 1 and B = Q was studied 
by |Shu| ( |1977[ ), and a generalization oi B for an isother- 
mal gas was pursued by Whitworth & Summers ( 1985 1 for 



solutions with weak discontirmities (Lou & Shen 2004). 

Clearly, for all n > and B / 0, Ba;'""^)/" is the 
leading term in the asymptotic solution of v{x) for large 
X. Specific for isothermal cases (i.e. n = 1 and 7 = 1), _B 
actually represents the constant reduced speed at a; -l-oo 
(e.g. Lou fc Shen|2004l. In reference to our previous work 



(e.g. Lou & Wang|2006 1, we shall refer to solutions according 



to their large— a; asymptotic envelope behaviours as follows 
in this paper: B = and t; > for "breeze solutions"; B — 
and V < for "contraction solutions" ; B > (thus v > for 
large enough x) for "outflow solutions"; B < (thus t; < 
for large enough x) for "inflow solutions" . 



^ In |Suto fc Silk| | |1988^ the term proportional to has been 
missed. Lou & Shi (2011 in preparation) first pointed out this 
missing term in the asymptotic expansion for very large x. 
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Similarity transformation ([5| and hydrodynamic PDEs 
lead to the following algebraic relation 



m — ax (nx 



(8) 



which implies the presence of the so-called "zero-mass line" 
(ZML) specified by na; — w = in the —v(x) versus x fig- 
ure presentation; below this ZML a negative enclosed mass 
M{r, t) would be physically unacceptable. 



2.2 The sonic critical curve and eigensolutions 

In two coupled nonlinear hydrodynamic ODEs ([gJi, we en- 
counter a sonic critical point when the denominators on 



the right-hand side (RHS) vanish ( Suto fc Silk 1988 n.b. 
n -I- 7 = 2), i.e. 



w = (2 - n) 



1/2 (1- 



1)12 



(9) 



In the three-dimensional (3D) variable space of x, v and a, 
condition ([9| represents a two-dimensional (2D) surface of 
singularity for two coupled nonlinear ODEs ([6|. In general, 
a solution cannot go across this surface smoothly due to the 
singularity of ODEs except for a special type of solutions, 
viz., the eigensolutions. One important necessary condition 
for searching such eigensolutions is that the numerators on 
the RHS also vanish simultaneously, which is equivalent to 
(see also |Suto fc Silk|1988| ) 

X — v){nx — v) 



, (nx — v) 



(10) 



While requirement ( |10[ ) is derived to have the numerator 
of da/dx to vanish, it is straightforward to verify that the 
numerator of dv/dx in eq. ^ also vanishes with eq. (|9| 
given (see |Suto fc Silk||1988) Lou & Shi 2011 in prepara- 
tion). In other words, conditions (|9| and (10 1 are sufficient 
to determine those eigensolutions crossing the sonic criti- 
cal curve (sec) smoothly, leaving one degree of freedom 
(DOF), i.e. the specific point at which the eigensolutions go 
across. ODEs (|9| and ( 10 1 are combined to determine the 
sec in the 3D variable space of x, v and a. The points at 
which the eigensolutions cross the critical surface smoothly 
are all located along the SCC. In the following, we shall use 
2:0, vq and ao to denote values of x, v and a on the SCC 
for solutions crossing the SCC smoothly. When solving two 
ODEs ([9| and ( 10 \ explicitly for xq, vq and ao, we find that 
the SCC consists of two segments represented by subscripts 
-|- and — respectively as two connected quadratic solutions: 



2:0, ± = 



/r, \l/2 {l-n)/2 

uo,± = na;o,± — (2 — n) ' Qq 
where 

Ka = n{n - 1) , 



(11) 



K,^i2- nf/' [n-l + ^) 4^"'/^ , (12) 
= 2(n - 2)aJ-" 



| |Lou fc Cao|2008| ). These two + and — segments of the SCC 
meet right at the point where 

ao = -3nV5n-2 + 2[2(4n-16nV21n^-9n'')l^''^ (13) 



We combine these two SCC branches to form a smooth 
global SCC. In the following, we refer to the segment with 
"-f" in eq. (Ill as segment 1, while that with "— " in eq. 
(Ill as segment 2, respectively. 

We then expand an eigensolution to the first order near 
the SCC as 



X = Xq -\- 5 



a = ao + a'(5 , 



Wo 



+ v'& 



(14) 



where a' and v' are the values of da/dx and dv/dx on the 
SCC for eigensolutions, respectively. Solving two ODEs ([9| 
and ( 10 1, we obtain a quadratic equation of v' as well as the 
equation of a' (see also Suto fc Silk|l988 1 : 

^(n+l)/2 



av' +bv' +d — , 



, xo — vo 



Xo 



(15) 



(2-71)1/2 

where three coefficients a, b and c are explicitly defined by 
a = 3 — n , 

ao 



b = 4{l-n)— +5n^ 

Xo 

2 



+ 



3n- 2 

(»-2) 
(3n - 2) 



2n- 



ao 



Xo 

2n + 4: . 



(16) 



Quadratic equation ( 15 1 for v' generally allows two roots. 



corresponding to two types of possible eigensolutions across 
the SCC at the same point on the SCC (iLou fc Cao|2008') . To 



be specific, with the formula of roots of quadratic equation 
( 15 1, we obtain the two roots of v' explicitly as 



-b ± {b^ - Aad) 
2a 



1/2 



(17) 



For an eigensolution crossing the SCC only once, those with 
v'_ will be referred to as type-1 solutions while those with v"^ 
are type-2 solutions. In general, those sonic critical points at 
which the solution take v' — v'_ will be referred to as type-1 
critical point, while v' = v\. for type-2 critical point. 

2.3 Shock jump conditions across the SCC 

For astrophysical applications, it is common to expect the 
existence of shocks in the hydrodynamic evolution of a gas 
sphere. Across a shock front propagating in a conventional 
polytropic gas, conditions of the mass conservation 



(P2 — Pl)Us = U2P2 — Ulpl , 

the radial momentum conservation 

Us{u2p2 — Uipi) = — pi +P2 ~ ulpl + ulp2 

and the energy conservation 



Pi 



P2 



■ U2 



7 — 1 7 — 1 

IP2 
7-1 



ulp2 



ulpi U2P2 
2 2 

IPi 



Ul 



7-1 



+ 



u'ipi 



(18) 



(19) 



(20) 



need to be satisfied in the co-moving shock reference frame- 
work in order to attain a physically acceptable shock discon- 
tinuity. We note here, especially for the energy conservation 
across the shock front, that the polytropic index 7 should re- 
main the same across a shock front for similarity solutions; 
if those 7 are different, n would be also different because 
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of n + 7 = 2. This would not be acceptable for a global 
similarity evolution in which a shock evolves self-similarly. 

It is useful to point out that those three shock equa- 
tions above are invariant under a commutation of subscripts 
1 and 2. We need to know the direction for entropy increase 
for identifying physical shock solutions. To be specific, we 
denote subscript 1 for the upstream flow. Solving combined 



shock conditions (18l — (20l and adopting self-similar trans- 
formation ([5|, those conservation conditions become explic- 
itly (a;si is the self-similar location of the shock on the up- 
stream side and Xa2 being that of the downstream side; in 
fact, they correspond to the same shock radius rs{t) but 
with different k values for the sound parameter) 



a2 



"1(7 + 1) (^1 - nxsi) 



X32 = Xsl 



"2(7 + 1) 



2(v — nXsi)^ — 
2nxsi + ui(7 — 1) + 



"7(7 - 



1/2 



(21) 



V2 = 



(^1 



■ Xsl Xs2 



(7+I) 



Xsl 



which are sufficient to determine those important self-similar 
variables {x, v, a\ at the immediate downstream side of 
a shock front when those of the immediate upstream side 
are known. The other way around, it is also straightforward 
to determine physical variables of the upstream side when 
physical variables of the downstream side are specified. 

In addition, we introduce parameter rjk for convenience 
as the ratio of different sound parameter k values at the 
upstream to downstream sides across a shock front, namely 




Figure 1. Examples of conventional polytropic self-similar so- 
lutions | |Lou fc Wang|p006[ l. These solutions have n = 0.7 (thus 
7 = 1.3). For a clearer illustration, the horizontal x-axis is shown 
in a logarithmic scale. The dash-dotted curve is the SCC, shown as 
a reference here. The light solid curve marked by "Lou & Wang" 
is the solution that crosses the SCC twice. It has two critical 
points on the SCC at (a, x, v) = (164.69, 0.002978, -2.4497) 
for a type-2 crossing and (a, x, v) = (0.1639, 3.4684, 1.5586) 
for a type-1 crossing, respectively (see |Lou fc Wang| ( |2006^ for 
more details). The heavy dashed curve marked by "Contraction 
- Free Fall" is the solution with contraction envelope and a free- 
fall collapsing core with asymptotic behaviour A = Ac + 0.3 (Ac 
is given by e g. l|28[ l), a polytropic counterpart of similar isother- 
mal results in |Shu| | [l977[ l. The heavy solid curve is the polytropic 
EWCS solution (A — >■ A'^). The dotted curve is the zero-mass 
line (ZML with nx = v). 



k2 



Xs2 
Xsl 



?(7 + l) 



1 



(22) 



2{v — nxsi)^ — aj{'y - 

Apparently, rjk refiects the strength of an expanding shock 
in a self-similar dynamic evolution. 

This rjk parameter here is equivalent to 1/A in the for- 
mulation of Lou & Wang ( 2006 1 . We now discuss the vari- 



ation of specific entropy along streamlines, reflected by the 
variation of k (or k; k = fc(47rG)'^~^i^'"^^~^' for a gen- 
eral polytropic gas and n = fc(47rG)'^~" for a conventional 
polytropic gas) across an outgoing shock front. The Mach 
number A4i of a shock is defined by 



Mi = \us — Ui\/si 



(23) 



where i — 1, 2 refer to immediate upstream and downstream 
sides, respectively, and Ui is the radial fiow speed near the 
shock front, Si is the sound speed, Us is the expansion speed 
of the shock front in our frame of reference. Using similarity 
transformation ([5|, we derive the Mach number of a shock 
in self-similar form on either side of a shock front as 



following discussion, we assume 7 > 1. It is then straight- 
forward to show that when Mi > 1, we must have z < 1 
Since z < 1 gives (see also Lou fc Wang|2006 1 

d(A^) 
dz 



d 

dz 



7(7'-l)(2-l)' 



[(7+1)^ -(7-1)]' 



< 



(26) 



and A — )■ 1 when 2: — )■ 1, this derivation yields that when 
Ml > 1, inequalities A > 1 and rjk < 1 (thus ki < k2 
and Ki < K2) will hold. This indicates the fact that the 
shock wave is supersonic in the upstream flow (A^i > 1) 
and causes an increase of entropy across a shock front (i.e., 
K2 > Ki). If however A^i is less than 1, subscript 1 would 
therefore denote downstream and the condition of entropy 
increase still holds (i.e., K2 < ki). In summary, shock jump 



conditions (I81 — (20l and naturally their solutions (21 1 are 



invariant with respect to a commutation of subscripts 1 and 
2; the irreversibility of a shock or the increase of specific 
entropy across a shock front should not be violated. 



M^ = 



7"7" 



where notations Ti and z are defined by 
= n - Vi/xi , z = r2/ri , 



(24) 



(25) 



which are the same as those in|Lou fc Wang| (|2006||.^ In the 



There is a typo in the definition of F, 
which is now corrected here. 



in |Lou fc Wang| | |2006| 



3 CONVENTIONAL POLYTROPIC 
SELF-SIMILAR VOID SOLUTIONS 

3.1 Several pertinent self-similar solutions 



Lou &j Cao ( 2008 1 



presented in Lou &; Wang ( 2006 ) and 



We first illustrate here several related self-similar solutions 
for conventional polytropic gas flows already obtained ear- 
lier. These are important in two aspects, viz., validating our 
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numerical code and preparing for further model analysis for 
self-similar voids surrounded by various dynamic envelopes. 

3.1.1 Conventional polytropic expansion-wave 
collapse solution (EWCS) 

Expansion-wave collapse solutions (EWCSs) were studied 
originally in an isothermal gas by Shu (19771. The natural 



generalization to conventional poljrtropic EoS was done by 
Cheng (19781. With a general polytropic EoS, this kind of 



solutions with 7 = 4/3 and n -I- 7 7^ 2 has been constructed 
by Lou & Cao (20081 in recent years. 



Conditions for EWCS can be readily derived from self- 
similar hydrodynamic ODE (|6|; meanwhile, those for the 
outer portion of a static envelope of a singular polytropic 
sphere (SPS) have been constructed (e.g. Lou & Wang 2006'). 
We can independently use some practical methods to obtain 
these solutions numerically with asymptotic solutions ^ in 
the regime of large x. As the coefficient of the leading term 
in the asymptotic expression of v{x), speed parameter B 
must vanish for this type of solution for a static polytropic 
envelope. For mass parameter A, we can get it from the 
asymptotic solution of v{x) by making the coefficient of the 
term of order a;^"'"^)/" to vanish, i.e. 



0" 



(27) 



2(2 -n) n ' 

uA" (3n- 2)_ 

which then yields {Ae is the limit for mass parameter A to 
achieve the conventional polytropic EWCS) 

l/n 

A^At, A,^ '"T^ ^' . (28) 



For the isothermal EWCS of Shu (1977), we have n = 1 and 
A^ = 2. This Aa value of A represents a limit. An exact 
equality oi A = Aa would correspond to the SPS solution. 

EWCS can be regarded as a special case bounding "col- 
lapse solutions without critical points" with A > A^ (e.g. 
Shu 1977 Lou & Zhai 2009[ ). Here we present one of such 
solutions with n — 0.7 in Fig. [l] (the conventional polytropic 
relation n -I- 7 = 2 holds). We found that the behaviours 
of solutions are qualitatively similar to isothermal cases as 
A gradually goes down to = 2 (see |Lou fc Cao||2008[ ); 
here the limit of A becomes A^ — 0.4044 for n = 0.7. For 
an exact conventional polytropic EWCS, a discontinuity in 
the first derivative appears at x = 1.3085, where the cen- 
tral infall side solution is tangent to the SCC. This discon- 
tinuity in first derivatives denotes the front of an expan- 
sion wave in self-similar evolution inside of which a collapse 
takes place and beyond which the gas remains a static con- 
ventional polytropic envelope (i.e. outer part of a SPS). We 
shall further construct this kind of solutions in the following. 

3.1.2 Solutions crossing the SCC smoothly 

We present here several eigensolutions crossing the SCC 
smoothly. In our Fig. [l] the solution curve marked with 
"Lou & Wang" is such an eigensolution example taken from 
figure 1 of |Lou fc Wang] p006) . This solution of n = 0.7 
(thus 7 = 1.3) is the conventional polytropic counterpart of 
isothermal envelope expansion core collapse (EECC) solu- 
tions of |Shen fc Lou| (|2004). Specific parameters are sum- 
marized in the caption of Fig. [T] 
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Figure 2. Conventional polytropic void solutions crossing the 
SCC smoothly with n = 0.67 (thus 7 = 1.33). Panel A il- 
lustrates the profiles of negative reduced radial speed —v{x) 
versus x and Panel B shows the corresponding profiles of re- 
duced mass density a{x) versus x where the a-axis is in the 
logarithmic scale. The dash-dotted curves in both panels are 
the SCCs, while the light dotted line in Panel A is the ZML. 
Three solutions of such type, marked by Model SI to S3 (S 
for "smooth") are plotted in light solid curve, heavy dashed 
curve and heavy solid curve, respectively. Model SI crosses the 
SCC at (xo, VQ, ao) = (4.277, 2.254, 0.02157), Model S2 
at {XQ, vo, ao) = (5.010, 2.736, 0.02358) and Model S3 at 
{xo, vo, ao) = (6.084, 3.442, 0.02657). Specific data and pa- 
rameters for each solution are labelled in the figure. 



3.2 Dynamic void solutions and their 

behaviours near the void boundary 

Among our self-similar solutions, it is possible for some to 
reach a specific line nx = v which is referred to as the ZML. 
By eq. the ZML leads to 

m\nx=v = ax^ijix — u) = , (29) 



and thus, according to eq. ||5|, 

^3/2j3n-2 







(30) 



' (3n- 2)G ' 

In other words, if there is a point in a solution where v equals 
to nx, the enclosed mass M within the corresponding radius 
oi nx = V (this radius expands with time t) would vanish. 

On one hand, these would suggest that there is an ex- 
panding void, inside which mass and gravity could be ne- 
glected as compared with those of the surrounding gas shell, 
embedded inside a time-dependent and self-similarly evolv- 
ing interface characterized by nx — v. On the other hand, 
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Figure 3. The relation of x^^ (wiierc Ocd a-nd v^^ can be deter- 
mined) versus xq for tiie conventional polytropic case of n = 0.67 
and 7 = 1.33. Tlie light solid curve marked along with crosses 
("x") corresponds to a branch of type-1 solutions, which is the 
only possible solution branch for xq < 50. That is, there are no 
type-2 solutions to reach x^d for xq < 50. 



condition nx — v implies that the expanding interface rep- 
resents a contact discontinuity (e.g. Lou & Zhai 2009, 2010) 
as revealed by the following relation 



xt 



(31) 



dt 

at the time-dependent interface radius where nx = v. This 
is one of the necessary conditions for the existence of contact 
discontinuity (e.g. |Lou fc Hu 2010 \ and the other require- 
ment is the balance of pressures on both sides, which will be 
discussed in Section |4] presently. In the following, we invoke 
subscript "cd" to refer to those variables at the immediate 
outside of the contact discontinuity (for independent self- 
similar variable x as an example, a;cd is simply the similarity 
location of the contact discontinuity surface). For instance. 
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X 

Figure 4. Conventional polytropic void solutions crossing the 
sec smoothly with n = 0.9 (7 = 1.1). Panel A illustrates the pro- 
files of negative reduced radial speed —v{x) versus x and Panel B 
shows the corresponding profiles of reduced mass density a{x) ver- 
sus X. The dash-dotted curves in both panels are the SCC, while 
the light dotted line in Panel A represents the ZML. Two solutions 
of such type, marked by Model S4 and S5 (S for "smooth") are 
plotted in heavy solid curve and heavy dashed curve, respectively. 
Model S4 is a type-2 solution crossing the SCC at (xq, vq, ao) = 
(4.571, 3.093, 0.66497) and Model 85 is a type-1 solution crossing 
the SCC at {xq, vq, ao) = (4.036, 2.608, 0.67397). 



We note here several properties of solution adjacent to the 
void boundary from the gaseous envelope side. According to 
self-similar hydrodynamic ODE Q, we can determine the 
first derivatives at the ZML, viz. 



dv 
dx 



2(1 -n) , 



da 

da:: 



n(l ■ 



(2- 



-a^^Xcd 



(33) 



and therefore the Taylor series expansions near the ZML 
are (this is actually the case of g = for a conventional 
polytropic gas in |Hu fc Lou||2008| |Lou fc Hu|2010 l 



v{x) = nxcd + 2(1 — n){x — Xcd) + ■ 



a{x) = Qcd + 



n{l 



-a"dXcd{x — a^cd) + ■ 



(34) 



Apparently, da/dx at the contact discontinuity vanishes if 
«cd ~ 0. This in turn hints that a{x) might be zero ev- 
erywhere, and we have indeed verified this by extensive nu- 



merical explorations (see also Lou fc Hu 20101. Moreover, 
leads to a singularity for the second order deriva- 



Q!cd 



tive of v{x) there (see Appendix A). In short, cases with 



Qcd = are not physically acceptable for a conventional 
polytropic gas. Meanwhile for rricd ~ 0, we must have 



lim a — , 



(35) 



indicating the presence of a jump or "clifF' in mass density 
across a contact discontinuity (this is exactly what the term 
"contact discontinuity" implies). The difference in mass den- 
sities inevitably leads to diffusion of gas particles across 
the contact discontinuity surface, which we shall discuss in 
Appendix [D] We shall focus on void solution cases of non- 
vanishing Ocd > 0. 



3.3 Conventional polytropic void solutions 
crossing the SCC smoothly 

Conventional polytropic void solutions can cross the SCC 
smoothly without shocks. As shown in subsection |2.2[ we 
are able to construct this kind of eigensolutions with central 
voids in self-similar evolution of dynamic expansion. 
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Panel Al n=0.815 
Xp=0.1 



Figure 5. The relation of x^^ (wiiere Oj-d a-nd '^cd can be deter- 
mined) versus xq witli n = 0.9 (7 = 1.1). Tlie light solid curves 
marked along with crosses ( " X " ) represent three branches of type- 
1 solutions, while those marked along with circles ("o") repre- 
sent one branch of type-2 solutions. Type-1 solutions are located 
in these intervals of xq: 0.23 < a;o < 0.66, 1.56 < xq < 2.19, 
2.62 < xo < 5.25. Type-2 solutions, on the other hand, are lo- 
cated in one interval of xq, i.e. 4.46 < xq < 5.25. The two right- 
most branches of solutions in different type actually converge at 
xq — 5.25. In fact, it is the very point where the square root in 
eq. \17\ vanishes and thus the two types of solutions are identical. 



While usual schemes always involve plotting a. a — v 
phase diagram by integrate towards a chosen meeting point, 
this series of void solutions can also be constructed as fol- 
lows without resorting to the a — v phase diagram. First we 
choose a point {xo, «o, Q^o) as the sonic critical point on the 
sec which satisfies ODEs ^ and (10 1. In fact, this repre- 
sents a DOF. There, specific types (either type-1 or type-2, 
depending on which of the two v' we take) of v' and a' are 



calculated consequently by solving eq. (15 I. After choosing 
a proper 5 < 0, we calculate v and a a.t x = xo + 5 < xo and 
integrate inwards to see if it could reach the ZML; when this 
is fulfilled, let 5 > to obtain the reduced variables 11 and a 
Sit X = Xo + S > Xo, we then integrate outwards to obtain a 
global solution for the self-similar dynamic evolution of an 
expanding central void. 

In general, the SCC and ZML together enclose an area 
in the figure for presenting v versus x profiles. Obviously, it 
is necessary for a void solution to have part of itself inside 
this area since it should touch the ZML. Therefore, there is a 
necessary condition for a solution crossing the SCC smoothly 
to reach the ZML as 



V < 



dvo 
dxo 



(36) 



where v' carries the same meaning as in eq. (15 I, of a spe- 
cific type (either type-1 or type-2). If this inequality is not 
satisfied along the whole SCC, we would expect no void 
solution that crosses the SCC smoothly of this type. It is 



straightforward to derive from eqs. (Ill and (15 1 that the 
necessary condition for the existence of type-2 void solution 
(with smooth behaviour across the SCC) is n > 0.840 at 
segment 1 of SCC (for the definition of "segment 1" , see eq. 
( |11[ ) as well as related discussions and definitions). By ex- 
tensive numerical investigations, we find no type-2 solutions 
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Figure 6. Presentation for the detachment (i.e. no cross-over) 
of the two phase curves in the phase diagram when n decreases 
to 0.80 (for the possible presence of shock solutions with EWCS 
envelopes). This detachment does not occur even when n is only 
slightly greater than 0.80 (say, 0.803). To be specific, we compare 
n = 0.80 cases with n = 0.815 cases in this presentation. In 
all four panels, dashed curves marked with circles ("o") are the 
phase curves representing different Xa\. The solid curve marked 
with crosses ("x") are the phase curves representing different 
Ocd a-iid Xcd — > 0+: these curves are the inner envelope of the 
"phase nets" with various combinations of a^d and aj^d- Panels 
Al and Bl illustrate the change of the phase curves in the "upper 
left region" (mainly for x^i < Xe) while n decreases from 0.815 
(Panel Al) to 0.80 (Panel Bl). Panels A2 and B2, meanwhile, 
present this change in the "lower right region" (for x^i > Xc) 
in the phase diagram. The detachment of the two phase curves 
(hence the detachment of the circled phase curve from the inner 
envelope of the "phase net" ) is apparent from this presentation. 



with Xo less than 50 for n = 0.67, whereas for n = 0.9, both 
types of void solutions exist. 

As examples, we present several solutions (n = 0.67 and 
n = 0.9) crossing the SCC smoothly in Figs. [2] and |4] with 
the most important relevant parameter xo being chosen as 
4.277, 5.010 and 6.084 for Model SI, S2 and S3, respectively. 
All those with n = 0.67 are type-1 void solutions, while there 
are both types of void solutions for n = 0.9. Corresponding 
diagrams for the relation of Xcd versus xo for n — 0.67 and 
n = 0.9 are displayed in Figs. |3] and [5] respectively. 

3.4 Expanding void solutions with shocks 

Shocks are another way by which solutions can go across 
the sonic critical surface. There is an extra DOF for void 
solutions with shocks: the shock location in the self-similar 
expansion. We have applied two types of numerical schemes 
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to construct void solutions with shocks. The following are 
the general outlines of these two procedures, perhaps with 
minor modifications in dealing with specific integrations. 

One computation procedure is to integrate inwards. For 
solutions with vanishing |w| and a at a; — >■ oo under consid- 
eration, we can use asymptotic solution ([7| to determine the 
initial value of integration at an x that is large enough (say, 
X = 30). We also choose a self-similar shock location in the 
upstream flow at a;^!. Then we integrate inwards numerically 
using the standard fourth-order Runge-Kutta scheme, apply 



shock jump conditions (21 1 in self-similar form at a: = a;si, 
and continue to integrate inwards from Xa2 and see if the 
solution would gradually approach and eventually touch the 
ZML. When such a solution touches the ZML, a global void 
solution with shock can then be readily constructed. Mean- 
while, values of dimensionless variables at the contact dis- 
continuity, Xcd {vcd = nxcd) and Ucd, are determined in a 
consistent manner. 

The other computation procedure is to integrate out- 
wards. Choose the self-similar location of void boundary x^d 
and reduced mass density Qcd there, and decide the location 
Xb2 of shock on the downstream side. Then integration is 
taken outwards to a relatively large value of x with respect 
of the jump at the shock by applying shock conditions (21 1 
(note the commutation symmetry with regard to subscripts 
1 and 2 discussed in subsection 2.31. The mass and speed 



parameters A and B to characterize the asymptotic solution 
behaviours at large x can also be evaluated with the value 
of v{x) and a{x) at a very large x by using eq. ([7|. 

Practically, there is not much difi'erence between the 
two kinds of void solution construction procedures except 
when parameters chosen in highly "sensitive" regimes. We 
have applied both numerical integration procedures as cross- 
checks of each other for the reliability of our void solutions. 



3.4-1 Expanding void solutions with EWCS envelope 

The outer part of EWCS is identical to the static SPS solu- 
tion, namely (e.g. |Suto fc Silk|1988| with n + 7 = 2) 







2(2-n)(3n- 2) 



l/n 



-2/n 



(37) 



beyond the expansion wave radius {xa for the similarity lo- 
cation of this radius), i.e. x > x^. For n — 1, this type 
of polytropic solution becomes the outer part of a singular 
isothermal sphere (SIS; e.g. |Lou fc Zhai||2009| 2010). The 
asymptotic behaviour of this kind of EWCS at small x is 
given by (e.g. |Suto fc Silk|1988[ |Lou fc Wang|2006 l 



2m(0) 
{3n-2)x 



-1 1/2 



(3n - 2)m(0) 



-,1/2 



2a;3 



(38) 



where constant m(0) represents the central mass point. 

We have obtained shock solutions with central void and 
EWCS envelope. When a specific envelope is prescribed, 
there is only one DOF, viz., the shock location. The tech- 
nique of V versus a "phase net", a variation of v versus a 

"phase net" 



phase diagram scheme IHunterl 



was developed in Lou fc Zhai (2009 



1977|) and the 

2010)], is applicable 



because of the similar situation. Steps of the scheme are 
summarized below. We first obtain the EWCS envelope by 
numerical integration. Choose a meeting point xp and the 




Figure 7. Void solutions with shocks and EWCS envelopes; Panel 
A for the negative reduced radial velocity —v{x) versus x profiles 
and Panel B for the reduced mass density a{x) versus x profiles. 
These conventional polytropic solutions have n = 0.9 (7 = 1.1). 
The SCCs in both panels are presented by dash-dotted curves, 
while the ZML in Panel A is shown by dotted line. Models El ("E" 
for "EWCS"), E2 and E3 are shown by light solid curve, heavy 
solid curve and heavy dashed curve, respectively. The envelopes 
beyond the shock of Models E2 and E3 are the same: SPS (thus 
Xsi > Xe) solution, but the locations of the shocks are different. 
Model El, on the other hand, has Xbi < Xe- 



similarity upstream location of the shock a;si. With shock 
jump conditions applied for connecting Xsi and Xs2, we nu- 
merically integrate from Xs2 to xf and record the value of 
V and a there. Select a series of Xsi values, get a series of 
{v, a) pairs at xp and use these {v, a) pairs to plot a curve, 
on which different points correspond to different values of 
3;si, in the phase diagram. Then we select a series of a;cd 
and Qcd in pairs and do similar things: integrate towards xp 
and get a sequence of curves or a "net" [i.e., the so-called 



phase net" in Lou fc Zhai ( 2009 1] in the phase diagram 



corresponding to different Xcd and Qcd in pairs. 

Applying this scheme, we have explored various void so- 
lutions with static polytropic envelopes. We find numerically 
that when the scaling index n is less than 0.80, there appears 
to be no such type of solutions with one shock. We realize 
for any specific n that the curve representing Xsi is actually 
bounded in the phase diagram and "shrinks" as n increases. 
Meanwhile, the "net" with different a;cd and Ocd pairs has 
an inner envelope (actually this envelope is outlined by the 
curve with different Ocd and Xcd — >■ 0^). When n becomes 
greater than 0.80, there are two regions in the diagram where 
the a::ai-curve intersects the "net" : one corresponds to the so- 



10 Yu-Qing Lou and Lile Wang 



lutions with Xsi < Xe, and the other, Xsi > Xe (in the phase 
diagram, these two regions are clearly separated). However, 
when n becomes less than 0.80, the Xsi-curve detaches from 
the net and neither of the two intersecting areas can still 
exist. We have tried various values of xf from 0.01 to 20 
and only find the same threshold of n — 0.80. 

We shall refer to the region with x^^ < x^ as the "upper 
region" and that with Xsi > Xe as the "lower region" in the 
caption of Fig.[6]because in a more complete phase diagram, 
the "upper region" always appears on the upper-left corner 
while the "lower region" is always on the lower-right corner. 
These terminologies can be readily clarified by referring to 
figure 7 of |Lou fc Zhai| ( |2009[ ) for an isothermal case. In 
addition, if there is no EWCS-envelope void solution (with 
shock), there will be no breeze solution with a void and a 
shock either. We have found that if < A < and B = 
(hence u > at large a:), the curve for Xsi will "shrink" 
even more, further preventing Xsi curve from intersecting 
the "net". Furthermore, the decrease of B parameter also 
leads to a "shrink" of a^si curve, therefore contraction and 
infiow solutions become much more difficult to exist. 

This detachment in the phase diagram related to the 
scaling index n — 0.8 is illustrated by examples in Fig. [6] 
Note that this separation by n < 0.8 is valid for conventional 
polytropic cases only; there is no such separation for general 



polytropic cases as shown in Hu & Lou ( 2008 1. Note also that 



void solutions with EWCS envelope and shock may still exist 
if the solutions cross the sonic critical surface for more than 
once (crossing smoothly or by shock discontinuity). 

Several void solution examples of this type are presented 
in Fig. [t] with n = 0.9 (thus 7 = 1.1) and Xsi being 0.4, 2, 
and 2.5, respectively. 



3.4-2 Void solutions with various dynamic envelopes: 
breeze, contraction, outflow and inflow 

With relevant parameters summarized in Table [ij a series 
of self-similar void solutions are obtained^ with shocks and 



n = 0.67 (values of n are to be discussed in subsection 4.2 1. 

Models 1 through 4 are four kinds of void solutions: 
those that have breeze, contraction, outflow and infiow en- 
velopes, respectively, all with n — 0.90 > 0.80. Their so- 
lutions are illustrated in Fig. [S] to show in non-dimensional 
form the profiles of reduced radial velocity v{x) and reduced 
mass density a{x) (definitions of these four kinds of en- 
velopes are given in subsection 2.11. Models 5 through 7 



are void solutions with outflow envelope beyond the shock 
with n = 0.67 < 0.80, which we shall elaborate for further 
applications. These solutions are displayed in Fig. [9] 

In numerical explorations for cases of n = 0.67, we have 
found that even for strong outflows, there are impressive 
trends to fall inwards at some intervals of x: this is clearly 
the results of self-gravity. This point is unambiguously illus- 
trated by Models 5 through 7, and other models not shown 



Table 1. Parameters specifying conventional polytropic void so- 
lutions with shocks. The first four models all have n = 0.9 while 
the last three all have n = 0.67. These solutions are constructed 
by integrating inwards and we confirmed their correctness by in- 
tegrating outwards within the tolerance of numerical errors. 



Model 


n 


A 


B 




Ilk 




"cd 


1 


0.9 


3 


-1 


1.5 


0.826 


0.809 


13.9 


2 


0.9 


3 


1 


3.5 


0.902 


2.509 


1.01 


3 


0.9 


2.542 





2.5 


0.905 


1.559 


2.03 


4 


0.9 


1.542 





2.5 


0.950 


1.204 


0.88 


5 


0.67 


2.43 


13.6 


5.5 


0.214 


0.239 


59.6 


6 


0.67 


2.25 


13.2 


7 


0.300 


1.031 


12.6 


7 


0.67 


4.50 


16.7 


9 


0.145 


0.358 


28.2 



here also appear similarly in this aspect. We shall further 
discuss this feature in subsection 14.41 



4 ASTROPHYSICAL APPLICATIONS 

There are several astrophysical situations where hot tenuous 
bubbles exist, shaped up by pertinent physical mechanisms. 
Here, we discuss cases that might be capable of describing 
the evolution of supernova at an early stage (the so-called 
"optically thick" stage), during which the predominant driv- 
ing power inside the stellar envelope could possibly be the 
pressure of some extremely relativistic particles such as neu- 
trinos and photons trapped inside a central cavity. Some 
necessary aspects need to be addressed before outlining our 
model scenario. 



4.1 Formation of a central bubble or cavity 

Research works have hinted at a scenario that a cavity can 
be formed surrounding the centre of a supernova during the 



initial phase. As concluded by Bethe ( 1990|, accelerating in- 
fall of substances is inhibited by powerful neutrino pressure 
from the neutrinosphere before neutrinos decouple from the 
gas and escape; moreover, a rebound shock may be revived 
when the intense neutrino flux is absorbed by surrounding 
nuclear matters at a radius r ~ 100 km or more. Such a 
mechanism can drive materials around the core outwards 
and shape up a rarifled bubble or cavity, which is clearly 
iUustrated in figure 2 of [Bethe fc Wilson] ( [T985| . This per- 
spective was further strengthened by simulations of |Janka"fcl 
pJillebrandt, ( ,1989a^b) and Janka & Miiller (1996), revealing 
the formation of a bubble around the centre of a supernova, 
filled with intense electromagnetic radiation field and other 
relativistic materials. 



^ The results in Table 1 are obtained with the inclusion of — (n — 
l)B'^x^~^/'^ /n term (Lou &; Shi 2011 in preparation). The values 
of mass and velocity parameters A and B, especially those with 
n = 0.67, have already been modified as compared with those 
without including this term. In contrast, n = 0.9 cases are almost 
not infiuenced by the inclusion of this term. 



4.2 Pressure balance, EOS, and the central power 

While it is allowed to be discontinuous in mass densities 
and temperatures across the contact discontinuity surface, a 
pressure balance across this contact surface should be main- 
tained to fulfill the necessary mechanical requirement. Here, 
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Figure 8. Presentation of negative reduced radial flow velocity 
—v(x) versus x (Panel A) and reduced mass density a{x) versus x 
(Panel B) for shock void solutions with inflow (heavy solid curve), 
outflow (heavy dashed curve), contraction (light dashed curve) 
and breeze (light solid curve) envelopes specified by Models 1 
through 4 in Table [T] respectively. The SCCs in both Panels A 
and B are shown by dash-dotted curves; the ZML is represented 
by the light dotted curve. 
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Figure 9. Presentation of negative reduced radial flow velocity 
—v{x) versus x (Panel A) and reduced mass density a{x) versus x 
(Panel B) for outflow-envelope solutions specified by Models 5—7 
in Table [l] The dash-dotted curves in both panels are the SCCs. 
The light dotted line in Panel A is the ZML. Among the three 
solutions. Models 5, 6 and 7 are illustrated with heavy solid curve, 
heavy dashed curve and light solid curve, respectively. The small 
inset at the upper-right corner of Panel A is a zoom of solution 
curves near the origin x = and —v = 0. 



we discuss consequences for the conventional polytropic EoS 
and some relevant aspects of the required pressure balance. 

For an extremely relativistic degenerate or hot gas, sta- 
tistical mechanics gives an adiabatic EoS (e.g. Callen|1985 1, 



(39) 



pp — constant . 

With a spherical volume V = 47rr'^/3, p is the mass density, 
p is the pressure, and the energy density is simply gotten by 
a multiplication of with c being the speed of light. 

Consider matters inside a highly rarified central cavity. 
For a temperature as high as ~ 1 — 10 MeV/fce or more, 
it is natural to expect substances (i.e. radiation field and 
products of electron-positron pair production processes) in- 
side the cavity to be relativistic. During the self-similar dy- 
namic expansion of a central void, transformation ([sjl gives 
Ted oc . Assuming that the mass (hence energy) inside 
the void is to be homogeneous and conserved (adiabatic) 
during the self-similar expansion, thermal pressure p is thus 



proportional to V 



-4/3 



(for relativistic matters, the 



speed is almost the speed of light c; hence homogeneity 
would be a sensible approximation for a not very large scale, 
say < c X Is). From those, we conclude that p oc f"*" for 
relativistic matters inside a void at the boundary (i.e., the 
surface of contact discontinuity). 



On the other hand in terms of the surrounding gas en- 
velope in self-similar evolution, the thermal pressure just 
outside the contact surface is proportional to t'^"~'^ by trans- 
formation ([5|; a contact surface has constant a^d and x^d 
in a self-similar dynamic evolution. There should be a sus- 
tained pressure balance across the contact surface, requiring 
at least the same time-dependence of pressures on both sides 
of a contact discontinuity. With this consideration alone, we 
would require n — 2/3 from 2n — 4 = 4n (i.e. 7 — 4/3 from 
n + -y = 2). 

We shall not consider the exact case of n = 2/3 for 
a homologous flow as done by |Goldreich & Weber] (|1980|) , 



Yahil (19831, and Lou & Cao (20081. Instead, we invoke a 



physically more plausible condition n = 2/3 -I- e with e > 
being a real number slightly greater than zero (e.g., n — 
0.67) to describe a small deviation from an adiabatic process 
for extremely relativistic matters inside a cavity, with a tacit 
assumption that the central compact remnant (e.g. a nascent 
neutron star or a nascent stellar mass black hole) continues 
to input energy into the "void" with a time-dependent rate, 
usually a decreasing one. 

Simple thermodynamic and mechanical considerations 
yield that an e added to n reflects some physical mecha- 
nisms adding more energy into a void during its dynamic 
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expansion, viz. 

dQ ~ 3epdV , (40) 

where Q is the amount of heat input into a void by some 
mechanisms as shown in Appendix B. 

By transformation ([5|, we readily derive 

^5/2^5e-5/3 



dQ 
dt 



3e 



+ e 



4/3- 



G 



(41) 



where G is the gravitational constant and k is the same 
sound parameter in transformation (|5|, and the subscript 
"cd" denotes quantities related to the contact discontinuity. 
The sound parameter k should be evaluated with in terms 
of dimensional values such as pcd, which depends on specific 
applications of our void model. 

We do not yet know the specific form(s) of energy re- 
lease from the central compact object into the surrounding 
cavity. We consider simple cases as examples of model con- 
sideration. 

There is a sensible physical possibility that the energy 
input comes mainly from the thermal radiation from the col- 
lapsing central compact object. As is conventional in study 
of neutron star cooling process, an "effective" temperature 
Tc is often introduced to indicate the ability of photon radi- 
ation (e.g. |Yakovlev fc Pethick|2004[ [Page et al.|2006l ). Note 
that even though neutrino cooling is much more important 
than photon cooling during this epoch, photon radiation is 
much more important for a sustained expansion of central 



void (see subsection 4.3 for details). With this Te specified, 
the radiation energy flux input into the bubble is given by 

^ ^ 47rr2aT^ 
at 

where rn is the radius of the proto-neutron star and a — 
5.6705 X 10"^ erg cm"^ s"^ is the Stefan-Boltzmann 



(42) 



constant. Combining eqs. (411 and (421 with n — 2/3 + e, 
we obtain a scaling of Tc as a function of t in the form of 

,5(n-l)/4 



Tc oc r 



(43) 



If simplified calculations of Page et al. ( 2006 1 hold even in 



the early evolution phase of a proto-neutron star, we obtain 
approximately for a "slow" ("standard") cooling process, 
Te oc i"^/^^ and thus n = 14/15 ~ 0.933, while for a "fast" 
cooling process, Te oc t~^^^ and thus n = 0.9. 

As an example, let us approximately calculate the 
amount of energy input rate with respect to the model ap- 
plied in our scenario for SN1993J in subsection |4.5.2| The di- 
mensionless model applied there has Xcd = 2 and Qcd ~ 
with n = 0.933 (in addition, the relevant Xs2 is 3.3). The 
dimensional counterpart is taken as what is specified in 
Fig. [TT] at t = 11 s. We then have the sound parameter 
k = 1.584 X 10^* c.g.s. unit and thus from eq. ([41 1 



^ ~ 8 X 10*^ erg s'^ 
dt ^ 



(44) 



A thermal luminosity as high as this would be possible in a 
very early stage in terms of pertinent calculations of |Page| 
et al. ( 20061) in their figure 14. For rn — 10 km, we obtain 



from eq. (|42[) an effective temperature Tc ~ 1.8 x 10^° K at 
that time, corresponding to a very early epoch of a proto- 
neutron star (see again figure 14 in Page et al. 2006 for a 
high surface temperature of a strange star). This tempera- 
ture Tc is much higher than the radiation temperature in the 



bubble (see Fig. 111. Moreover, the temperature in the bub- 
ble drops at a rate proportional to — t""'^^'*, which 
is much faster than Te oc t~^^^'^. This can be seen from the 
pressure confining the central bubble p = Pcd oc t^"-* [see 
also transformat ion ([s] )] and p on T"^ for radiation [see eq. 
(3.53) in |Callen| ( p85| ]. 

To continue the above consideration, we may presume 
approximately a black body emission spectrum from an 
isothermal proto-neutron star. If the heat capacity at con- 
stant volume of a proto-neutron star depends on the tem- 
perature in a power-law form, viz. Cv oc T* with ^ being an 
exponent, we have 



T^dTn oc d(7n oc -T^dt 



(45) 



where U-n is the internal energy of the proto-neutron star. 
From that, we immediately derive the proportional relation 
Tn oc or 



dQ 
dt 



df/n 

' dt 



oc Tt oc 



(46) 



In reference to eq. (41 1, we need to require 4/(^ — 3) = 
5e — 5/3, leading to the following relation between n and ^ 



n = 1 -f 4/[5(e - 3)] 



(47) 



Since for a degenerate Fermi gas, theoretical computations 
of heat capacity at constant volume Gv gives the exponent 
^ range of < ^ < 1 (e.g. Greiner et al.|[l995 from its fig- 
ure 14.2). It follows that scaling index n would be restricted 
to the range 3/5 < n < 11/15 according to eq. (471. Ac- 
tually, we should require n > 2/3 in our formulation for 
physical similarity solutions. In other words, the mechanism 
attributing the energy input to the black body radiation 
from a central proto-neutron star to a surrounding hot bub- 
ble or cavity can only be responsible for the cases where 
inequality 2/3 < n < 11/15 is satisfied. We emphasize that 
if for some reason, ^ is allowed to be negative, then the value 
of n may approach 1. For other physical mechanisms of en- 
ergy input into the central cavity from a proto-neutron star, 
inequality n < 11/15 may not be necessary. For example, 
numerical simulations of [Thompson et al.| ( |2001[ ) have of- 
fered a neutrino luminosity proportional to t~"'^ (see their 
figure 9) after a SN explosion and this would correspond to 
n = 0.82 > 11/15. 



4.3 Momentum transfer by scattering processes 

The initial rapid detachment of an out-fiowing envelope from 
a collapsing compact core is primarily accomplished by a 



powerful neutrino pressure (e.g. Bethe 1990). After hundreds 



of milliseconds, energetic neutrinos decouple from surround- 
ing gas materials and escape, an important source of driving 
power for the bubble expansion vanishes. 

Even though the radiation (mainly involving trapped 
photons and electron-positron pairs produced by pair pro- 
duction) is relatively weaker in luminosity than neutrinos, 
scattering process, by which the momentum and energy are 
transferred, of photons by matters in the gas shell is never- 
theless much more effective. 

Meanwhile, previous studies have implied that the 
amount of energy needed to blow a massive stellar envelope 
up is roughly at the same magnitude (~ 10*^ erg) as that 



radiated by photons in SN explosion (see Janka & Miiller 
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|1996| for details; this estimate could also be derived from an 
integration of the energy spectrum in Chevalier|1974 1. These 
also suggest that trapped photon radiation may be a major 
driving power of further detachment between a collapsing 
core and an expanding massive envelope after neutrinos have 
already shaped an initial central cavity and decoupled from 
surrounding gas materials. 

Such a model consideration is outlined in the following 
subsection, where we compare the contributions of neutrino 
flux and of photon radiation field. We will examine the possi- 
bility of photon radiation field of being an important driving 
power to the expansion of a shocked massive envelope and 
a sustained central void expansion. 

4-3.1 Scattering of neutrinos by heavy nuclei 

Energetic neutrinos are important in the initiation and evo- 
lution of SN explosions (e.g. Bethe|1990 1. Scattering process 
of neutrinos with extremely high density matters is treated 
as coherent scattering, with a scattering cross section 
given by (e.g. |Tubbs fc Schramm|1975| ) 



(48) 



where A is the number of nucleons in one nucleus, E,, is 
the neutrino energy, rUe is the electron mass and c is the 
speed of light. The electron rest mass energy nicC^ comes 
basically from the process n — >■ p + + i/c in a nuclear 
reaction. Equation (48 1 does not include ^ neutrinos i/^ (i?^) 



and T neutrinos f-r {i't), because the latter two are much less 
significant as compared with electron neutrinos Ve (fc'c). 

For a highly degenerate dense stellar core, the scattering 
of neutrinos are so effective that a central "neutrino bubble" 
on a spatial scale of ~ 100 km is shaped up by the extremely 
intense neutrino flux released from the collapsed core (e.g. 
Bethe|1990[|Padmanabhan|2001[ ). In this case, the mean free 
path of electron neutrinos is given by (e.g. Padmanabhan 
200l| ) 



3 X lO"* cm ( — 



10^^ g cm ^ 



-5/3 



(49) 



where A is the average number of nucleons per nucleus. On 
the other hand, the shaping of the neutrino bubble causes 
the matters surrounding the core to expand significantly, 
making them less dense and non-degenerate gradually (see 
Appendix |D]| . For a non-degenerate gas envelope, is es- 
timated by 



A. 



Amp 
E,.p 



0.5RqA 



E, 



1 MeV 



IQi" g cm- 



(50) 



This Xv would be too long to accomplish an effective mo- 
mentum transfer for a mass density less than that of nuclear 
matters. Therefore we expect the material shell to be trans- 
parent for neutrinos shortly after a "neutrino bubble" being 
shaped up. 



4-3.2 Scattering of photons by charged particles 

Classical and semi-classical theories yield almost the same 
result for the scattering cross section Eph between a photon 



and a charged particle, usually referred to as the "Thomson 



cross section" given by (e.g. Jackson||l999 1 



■Jph 



Stt 



(51) 



where e — 4.803 x 10~^° e.s.u. is the electron charge, m is 
the charged particle mass and c is the speed of light. 

Clearly, Thomson cross section is much larger for elec- 
trons than for nuclei because me <^ rUp (specifically for elec- 
trons, Eph — 6.65 X 10~^^ cm^). For a fully ionized plasma 
of mass density p, we obtain the mean free path of pho- 
tons (rinuc and Uc are the number densities of nuclei and 
electrons, respectively): 

ph,c Eph,G^c 

Numerically, we estimate a photon mean free path as 



A 



ph 



5 cm 



1 g cm^ 



(53) 



This implies that photons are tightly trapped inside the 
cavity even for an envelope mass density as low as ~ 
10~^ g cm~^. The much higher efficiency of momentum 
transfer (compared with neutrinos) also permits the pos- 
sibility of radiation-driven envelope expansion of a void. 



4.4 Void solutions as asymptotic conditions 

The mass is conserved in evolution and thus mass densi- 
ties throughout the envelope decrease as time goes on for an 
overall expansion. Since we investigate the situation that ra- 
diation field continues to drive the envelope outwards, when 
the gas in the inner envelope (especially near the contact 
discontinuity surface) becomes sufficiently tenuous, trapped 
photons can leak out gradually. 

We define At as the time duration after ti (for its defini- 
tion see Appendix C; hereafter we use subscript "i" for those 
initial values in dynamic void evolution) to the time that the 
mass density near the contact discontinuity beco mes fairly 
low (e.j 



4.3.21 



pcd ~ 10~ g cm~ , see subsection 
In reference to similarity transformation (|5l, an esti- 
mate of At is given by 



At ~ t. 



Pcd,: 



10~5 g cm" 



1/2 



- 1 



(54) 



Meanwhile, the diffusion across the contact discontinuity 
may modify our model. In fact, diffusion process is not so 
effective as analyzed in Appendix [D] 

The mass density near the contact discontinuity is at- 
tenuated in a self-similar expansion, which would make our 
radiation driving expansion model invalid after a sufficiently 
long lapse in time. We here show that our self-similar model 
may be further used as an asymptotic solution when the 
radiation-driving mechanism is no longer effective. 

Recalling self-similar transformation ([5|, the pressure 
decreases with pcd oc while the radius of the contact 

discontinuity expands with a power law i". We define the 
integral of radiation pressure pcd over the contact disconti- 
nuity surface as the "radial radiation force" Fr, 



ri /I 2 2 7 I 

Fr = 47rredPcd = a^cd^cd 



G 



(55) 
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Figure 10. Dimensional model presentation at t = ti (the initial 
timescale of a self-similar void evolution). Panels A, B and C 
illustrate the radial flow velocity (in unit of 10^ cm s~^), mass 
density (in unit of g cm~^) and temperature (in unit of Kelvin 
K) profiles of the model, respectively. Vertical line-sections in the 
plot on the left represents the contact surface, inside of which is 
the intense radiation field, while those on the right indicates the 
expanding shock. Note that we have also shown the temperature 
of the radiation field inside the central cavity in Panel C. 



Table 2. Parameters specifying conventional polytropic void 
models with shocks; in dimensional form, they describe the dy- 
namic evolution of a gas shell, e.g. a SN. "CD" in this Table 
denotes the contact discontinuity. "Duration time" is the time 
span of the model validity considering the scattering process of 
the charged particles by photons (see subsection |4.4| for details). 



Item 


Variable Name 


Value 


Total mass 


AI 


~ 20Mq 


Density at CD 


Pcd.i 


~ 4.7 X 10^" g cm 


Pressure at CD 


Pcd.i 


~ 10^^ dyn cm" 


Initial void radius 


^cd,i 


~ 160 km 


Radial velocity of CD 


^cd,i 


~ 6000 km s-l 


Cavity Temperature 




~ 9 X IQi" K 


Duration time 


At 


> 10^ s 



become more apparent, difTusion across the contact surface 
may become severe), void solutions are still reasonable to 
the degree of (at least) approximations, since the inertia be- 
comes predominant in further continuation of expansion. 

4.5 Self-similar evolution of SNe driven by 
central photon radiation pressure 

In this subsection, we briefly discuss a few specific applica- 
tions of the conventional polytropic self-similar void solu- 
tions with shocks in the context of SNe. The procedure of 
constructing a physical model from dimensionless solutions 
are summarized in Appendix C. 



This result indicates that the force exerted on the gas shell 
by the radiation inside the bubble becomes less and less im- 
portant cis time goes on, since there is usually 4(n — 1) < 
(from 7 > 1 and n -I- 7 = 2). For a long enough time t, 
the pressure force across the contact discontinuity surface 
diminishes to a negligible level. In other words, even though 
there should be a certain pressure across the contact dis- 
continuity surface, the dynamics of the shell are not sig- 
nificantly infiuenced if those pressures are weakened when t 
becomes sufficiently large. This analysis enables us to regard 
our void model to be valid as an asymptotic solution in the 
epoch when the gas shell becomes "optically thin" and thus 
radiation in the bubble leaks out. 

As an example, for the model parameters specified in 
Table [2j if we take pcd ~ 10~^ g cm""^ as the transition 
criterion from "optically thick" to "optically thin" , we get 

Pcd|p^^=io-5 g cm-3 = 3 X 10^ dyn cm"^ . (56) 

This is a low pressure in comparison with the inertia of the 
stellar envelope. We evaluate the "radial radiation force" 
exerted by the radiation pressure as 



Fr = 47rrcdPcd 



10^^ dyn 



(57) 



As a result, even though central radiation field has gone 
and there may be many deviations from the original forma- 
tion (e.g., the shock may eventually die out, asymmetry may 



4-5.1 A dynamic void model of self- similar evolution 

By properly specifying similarity transformation and enve- 
lope cutoff (see Appendix C), we would adopt our Model 6 
with key parameters summarized in Table [l] to describe a 
self-similar void evolution for a model SN. The progenitor 
has a mass of ~ 20Mq. An envelope cutoff radius is set at 
rcut ^ 40i?Q where the mass density is pcut — 10~* g cm~^ 
and the gas temperature is taken as T ~ 4000 K. 

Before the phase of self-similar evolution, a cavity 
around the center has already been carved out by the power- 
ful neutrinosphere, before neutrinos decouple from the sur- 
rounding envelope and escape. The cavity radius is initially 
estimated as > 100 km, by Bethe & Wilson (1985), Bethe 
(1990) and later works on the Wilson mechanism such as 
|Janka fc MMer| ( [T996| . Here we take the initial cavity ra- 
dius Tcd.i to be ~ 160 km. Subsequently, the stellar enve- 
lope surrounding the cavity expands outwards continuously 
powered by photon radiation pressure during the following 
self-similar evolution. Some important numerical values of 
parameters for describing this model (which also specifies 
the similarity transformation) are given in Table [2] 

While At (defined in subsection 4.41 is presented in Ta- 
ble [2j we note that the model can remain valid after this At 
as a asymptotic condition noted in subsection |4.4| 

As an example, radial flow velocity, density and tem- 
perature profiles at the initial stage {t — ti) are displayed 
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in Fig. |10[ From Panel A of Fig. |10[ we observe the impact 
of the radiation-driven mechanism in the dynamic evolution 
of envelope with an expanding central void. While the ra- 
dial flow speed throughout the shell remains always positive 
and large, there are still intervals of radius where gas ma- 
terials show an impressive trend to collapse inwards under 
self-gravity. In fact, a quick glance at some other self-similar 
models with 7 in EoS being close to 4/3 suggests a strong 
trend of the shell to decelerate (or even collapse) in the up- 
stream side of a shock (see Models 5 and 6). The expanding 
shock here serves as an important and effective accelerator 
to drive the envelope outwards against self-gravity. 

Regarding our analysis of self-similar model, the pho- 
ton radiation field trapped inside the central bubble appears 
to be a physically plausible candidate of the energy source, 
which drives the envelope outwards continuously during the 
optically thick phase. From Panels B and C of Fig. |10[ we 
readily recognize an outgoing shock. It is a strong shock 
with a drastic jump across the shock front in both mass 
density and temperature. The contact discontinuity surface 
at Tcd.i = 160 km initially, on the other hand, was supported 
by the strong neutrino flux before the decoupling of neutri- 
nos. The neutrino flux also heats up rarifled materials inside 
the void signiflcantly and is responsible for the discontinuity 
in temperature across the contact interface. 



4.5.2 Supernova SNl 993 J 
SN1993J in M81 is a type-IIL core collapse supernova about 



3 Mpc away from us (e.g. Schmidt et al.|[l993 ). This close 
neighborhood has enabled a plethora of observations. Here 
we take the observation and simulation results presented in 
Marti- Vidal et al.| ( [2011a|b[ ) and references therein. SN1993J 
has a spherical asymmetry < 2% and it would be of consid- 
erable interest to explore the capability of our void model 
with spherical symmetry. Numerical simulations of IMarti^ 
Vidal et al.| (|2011a|b[ ) have adopted the dynamic model of 
Chevalier ( 1982 1 where the self-gravity is neglected. 



There are different flttings of SN1993J expansion pre- 
sented in [Martf- Vidal et aT] ( |2011a[ ) with various expansion 
parameter m (not the reduced mass here) which takes the 
similar connotation as n in our formulation. An accepted 
selection of n seems to be n = 0.933 before ~ 360 days [sim- 



ilar to the selection specified in Marcaide et al. (2009)]. Here, 
t ~ 360 days is observed as the "break" time, after which an 
observable deceleration in expansion occurs (see subsection 



4.2 



for the implication of n). 
Some model features are noted here. Several simulations 
such as [Martf- Vidal et ah] ( |2011a[ ) have suggested that the 
radiation opacity of the shell has been fitted to be more 
than ~ 80% during 1 yr ~ 3 x 10^ s. This implies that our 
model may be valid during the first ~ 360 days. This justifies 
our void model to be applied even in later times since the 
pressure is actually low enough to be ignored dynamically 
(see subsection |4.4[ ). 

Here we briefiy describe a possible version of dimen- 
sional dynamic model for certain features of SN1993C ejecta. 
The mass of the progenitor is ~ 17M© (e.g. [Aldering et al.| 
). The radio bright shell consists of ~ 30% of the radius 
(e.g. Marcaide et al.|2009 1 (it is taken that this radio bright 
region consists mainly of the shocked gas between the shock 
and the contact discontinuity) and so forth. 
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Figure 11. Possible fittings of radial velocity {u in Panel A), 
mass density (p in a logarithmic scale in Panel B) and temper- 
ature (T in a logarithmic scale in Panel C) versus radius (r in 
a logarithmic scale) profiles of SN1993J ejecta in the early stage 
(a timescale of ^ 10 s). The "shell" region between the contact 
discontinuity and the shock is ~ 30% of the shock radius, in ac- 
cordance with radio observations (e.g. [Marcaide et al.|2069} . 



For our void model, the mass cutoff is set at r^u 



where the mass density is pcut ~ 0.03 g cm , according to 



the mass-radius relation in Padmanabhan (20011 



The "bright" region is assumed to be between the shock 
surface and the contact discontinuity where the gas is com- 
pressed and shock heated and charged particles (electrons 
in particular) are dramatically accelerated across the mag- 
netized shock front - this gives considerable radio emissions 
which can be detected by radio telescopes. Expansion of the 
shocked spherical shell is self-similar; specifically, the shock 
radius will expand by 



= 3.9 X 10' 



(58) 



where the timescale ti is fitted to be 11 s. If we input t to 
be 330 days after the SN explosion, we would expect the 
radius of the shock front (hence the radius of detectable 
radio emissions) to be = 3.8 x lO^'^ cm. Taking SN1993J 
to be 3 Mpc away and neglecting cosmological effects [e.g., 
[Ryden[p002| ], we would then have the angular radius 50 as 



3 Mpc 



0.85 mas 



(59) 



which is in very good agreement with the result of VLBI 
observations [~ (0.818 ± 0.015) mas] presented in [Marcaide[ 
et al. (20091. Fig. 11 shows these results for comparison. 
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5 SUMMARY AND CONCLUSIONS 

We have systematically explored various conventional poly- 
tropic (i.e. n + 7 = 2) self-similar void solutions by adopt- 
ing self-similarity transformation ||5| in nonlinear Euler hy- 
drodynamic PDEs of spherical symmetry. Different types 
of dynamic void solutions with outgoing shocks are con- 
structed numerically, while those with n — 2/3 + e and e > 
(e.g., n = 0.67) are described in more details. While highly 
idealized, these are dynamic void solutions with expanding 
shocks into various possible envelope types very close to the 
so-called "hot bubble" or cavity cases with potential astro- 
physical counterparts. Combined with proper central sources 
of energy and momentum, such void solutions can be initi- 
ated and sustained by matching the pressure balance condi- 
tion across the contact discontinuity in expansion during a 
certain phase of self-similar evolution. We discussed several 
possible situations. 

For supernova explosions of massive progenitors, we ad- 
vance the following physical scenario for the gross separa- 
tion between a more massive envelope and a less massive 
collapsing central compact object. Shortly after the onset 
of core collapse of a massive progenitor star, the powerful 
neutrinosphere trapped in an extremely dense core drives 
outwards to carve out a rarified zone of one or two hundred 
kilometers. As energetic neutrinos escape from the stellar en- 
velope of decreasing density, intense photon radiation field 
and/or electron-positron pair plasma mixed together con- 
tinue to drive the central cavity outwards against the en- 
velope. When the mass density of the expanding envelope 
becomes sufficiently low, photons leak out effectively in the 
optical thin regime eventually. If the pressure across the con- 
tact discontinuity interface diminishes to a sufficiently low 
level or diffusion effects (see Appendix D) smear out the 
"contact discontinuity" at sufficiently large radii, the outer 
remnant envelope continues to expand by inertia into the 
tenuous interstellar medium. 

Specifically for our model construction, we have first ex- 
plored possible conventional polytropic void solutions cross- 
ing the sec smoothly. Similar to those isothermal voids 
shown by |Lou fc Zhai| ( |2009[ 2010), there are several in- 
tervals of the independent similarity location of the sonic 
critical point xq — some of which permit the existence of 
such type of void solutions, while others do not. The SCC 
are divided into two segments for n ^ 1, and there are no 
such void solutions crossing SCCs at type-2 critical points 
in segment 1 of SCC for n < 0.84. In fact, there is no type-2 
void solution with xq < 50 for n less than 0.84, as we have 
found in extensive numerical explorations. 

Self-similar void solutions with outgoing shocks are also 
constructed, with several different kinds of envelopes outside 
the expanding shock, including the static SPS envelope, the 
outflow, inflow, breeze and contraction envelopes. Through 
extensive numerical explorations, we found no void solutions 
with shocks propagating into a static SPS, or breeze, or con- 
traction envelope for n less than 0.80. 

We have constructed some dynamic void solution ex- 
amples, especially in the context of SNe, as applications of 
our self-similar void models. For hot bubbles, 7 values are 
usually 4/3 (see Lou & Cao 2008 for a general polytropic 
gas) and for conventional polytropic gas dynamics, the cor- 
responding n should be 2/3 exactly. In reality, we expect 



deviations from the exact n = 2/3 cases. We have thus dis- 
cussed deviations of n values from 2/3, especially the phys- 
ical reality of e in n = 2/3 -I- e. In addition, we have dis- 
cussed the diffusion effect at void boundaries in Appendix 
D, flnding that the diffusion is around merely 1.3 percent as 
the central cavity radius nearly doubles. We have compared 
the roles of radiation fleld and neutrino flux respectively af- 
ter the neutrinosphere having decoupled from the gas shell 
with the result that radiation field would be a dominant 
force in the evolution of SN ejecta while the ejecta tend to 
be transparent to neutrinos. Also in subsection |4.4| we have 
discussed applications of our model as an asymptotic con- 
dition when the pressure across the contact discontinuity 
diminishes to insignificant level. 

Finally, we have attempted to specify our model param- 
eters to describe dynamic evolution of astrophysical objects 
and construct dimensional profiles in Fig. |10[ An example 
of application is the SN1993J. The dynamic model has been 
used to portray SN1993J evolution of expansion before the 
"break time" thr, where we have found that the results are 
able to capture some characteristics presented by radio ob- 
servations, such as a ~ 30% radius bright shell, the angular 
radius and the total stellar mass. 
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APPENDIX A: PROPERTIES FOR THE JUMP 
OF a FROM ZERO AT THE ZERO MASS LINE 

At the zero mass line (ZML) where nx — v, eq. Q shows 
a singularity for dv/dx and da/dx if we simply take a — 
there with 7 > 1. An alternative approach is to let a — > 0, 
which leads to 



da 
da; 







dv 
dx 



according to eq. ([6| 
to X can be calculated below as 

-]2 



->2(l-n), (Al) 

nx — v 

The second derivative of v with respect 



d'v 
dx^ 



(3n - 2) 



2 + 



n{n — l)x 
in- 2) ' 



(A2) 



It is readily seen that d^v/dx^ is equal to 2/x in an isother- 
mal gas with n = 7 = 1 (see eq. 23 of Lou fc Hu|2010 1 . Nev- 
ertheless with 7 > 1 (thus n < 1 in a conventional polytropic 
gas), we would have d^v/dx'^ — ^ 00 for a — >■ 0. Moreover, this 
a\nx^v — >■ would also give a solution with zero a every- 
where (as a function of x), which has also been verified by 



We can also attain these results by setting a = in related 



equations (16) and (17) of Lou & Hu (2010i 
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Lou & Hu (20101. Therefore, we cannot obtain a physically 
acceptable solution for a = at ZML in the conventional 
polytropic void model of self-similar evolution. 



APPENDIX B: THERMODYNAMIC 
DERIVATIONS OF ENERGY INPUT RATE 

In Appendix B here, all quantities are for the radiation field 
inside a bubble. From the energy conservation, we have 



dU ^dQ~ pdV , 



(Bl) 



where dQ is a differential heat transfer into a bubble. Ra- 
diation and relativistic matter have the internal energy 
U = 3pV where p is the pressure and V is the bubble vol- 
ume, and thus 

dQ = 4pdV + 3Vdp . (B2) 
If the EoS takes the form of p oc p*'^^"", we would have 



dp _ f 4 \ dp 



4 \ dV 



and it follows immediately that 
dQ = SepdV . 



(B3) 



(B4) 



APPENDIX C: FROM DIMENSIONLESS 
SOLUTIONS TO PHYSICAL MODELS 

In Appendix C here, we discuss how to produce a physical 
model from the reduced dimensionless solutions of ODE 
It is important to reasonably specify the self-similar trans- 
formation (there are some DOFs) and hence consistently in- 
terpret one dimensionless solution for a physical model. We 
can then determine the subsequent self-similar evolution of 
a void surrounded by a shocked envelope. 

A convenient means to consistently specify the transfor- 
mation parameters is to relate the initial values of physical 
quantities of the flow system on the contact discontinuity 
(e.g., the mass density pcd) to the corresponding dimension- 
less variables (e.g. Ocd) with respect to eq. ([5|. We shall 
attach a subscript "i" for those initial values of physical 
variables in a self-similar hydrodynamic evolution. 

First we consider the time t. With self-similar models 
adopted here, the time t cannot be zero even at the very 
first stage of similarity evolution, because t = is either 
a zero point (n ^ 1) or a singular point (0 < n < 1) of 
the self-similar transformation ([5| and resulting ODE 
A value ti, indicating the value of time t at the initial phase 
of self-similar evolution, needs to be specified or chosen. The 
self-similar transformation ([5|, from the initial mass density 
at the contact discontinuity surface pcd,i to acd, enables us 
to estimate ti as 



U = 



\47rGpcd,i 



1/2 



= 0.011 s a 



1/2 



Pcd.i 

IQiOg cm-3 



-1/2 



(CI) 



A similar ti is introduced as the "cut-off time" in |Lou fc| 



Wang ( 2006 \ with respect of shock evolution; alternatively. 



we here focus on various dynamic evolution of voids. 

We can also determine the value of sound parameter k 



at ti which is still valid in the subsequent evolution (for the 
presence of a shock in self-similar evolution, this k would 
then be k2 for the downstream flow) 



a/2 



— ± , - 



(C2) 



Another indispensable consideration deals with the cutoff 
enclosed mass and the cutoff mass density. For our self- 
similar void model, we need to choose a sensible enclosed 
mass M. For a very large r corresponding to a; ^ 1 at a 
fixed t, asymptotic solution a oc x~^^" (eq. ([7|) yields 



2/ N 3-2/n 

: ax (nx — v) nx ' 



(C3) 



As a; — >■ oo, m — >■ cxd for n > 2/3 and m — )■ for n < 2/3. 
Neither of them is physically acceptable as M is proportional 
to m and M/m is independent of radius in our self-similar 
solution. In dealing with this problem it is sensible to set 
a cutoff radius within which materials are included in the 
stellar mass; beyond this cut-off radius, materials are not 
regarded as stellar mass. In order to check whether the cutoff 
is proper, we can introduce the "cutoff density" - the mass 
density at the cutoff radius, as a criterion, which should 
be low enough for a reasonable cutoff. The cutoff radius is 
usually large but constant. 

With the help of asymptotic solution ([7|, we derive the 
following expressions for the enclosed mass at the cutoff ra- 
dius Tcut which is large enough: 



Mcut _ 2 ^ 10^^(1.44 X 10-^)^/" ^"^cd'^" ^l/n-3/2 



Ma 



(3n - 2)acd 



r-cut 



3-2/n 



IQfi cm 



2/n 



Pcd,i 

1010 g cm- 



and the cutoff mass density 



pcut 



109g cm 



- = 0.99 (1.44 X lQ-'' f-Anl^-x-J'-a-J 



rent 



-2/n 



/_rcd^\ 
V10« cm) 



f Pcd,, 



VlQiOg cm- 



(C4) 



(C5) 



Notations for A, n and k in these expressions carry the ex- 
actly same definitions as those in eq. ([sjl . Still, the sound pa- 
rameter k here is for the downstream side, fc2, in the presence 
of a shock. We demonstrate below that the enclosed mass 
at a large x is independent of time t. At a large radius r for 
a given time t, the enclosed mass is evaluated by 



,3n — 2 2/ \ 

* ax [nx — V) 



3n-2^^3-2/n 



OC f' 



-n\3 — 2/n ,0 
I ' CC t 



(C6) 



In the above derivation, we neglect v in (nx — v) because of 
a small v as compared with x in reference to eq. ([7|. The 
leading order term of x in eq. ([7| says, when x 2> 1, there 
is n ~ x^'^^"". Therefore v/x ~ and converges to 

zero when x ^ oo (for n > 0). This independence of time 
can also be explicitly seen from the fact that the radial flow 
velocity is usually small at large radii. Thus, the radial flow 
of materials at large r can be neglected in terms of its effect 
on the total mass enclosed by the spherical surface there. 
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APPENDIX D: DIFFUSION ACROSS THE 
CONTACT DISCONTINUITY INTERFACE 

The balance of pressures across a contact discontinuity is 
a necessary condition for the estabhshment of the disconti- 
nuity interface bounding an expanding void during the pro- 
cess of self-similar hydrodynamic evolution. Yet a dramatic 
difference of mass concentration (i.e. difference in chemical 
potential, to be precise) should lead to diffusion across the 
contact discontinuity interface. We now estimate diffusion 
effects across the contact discontinuity surface surrounding 
the central void - more in the context of a supernova ex- 
plosion. Somewhat different from what has been explored 



by Lou & Zhai (2009 2010), the envelopes or gas shells are 



at a high temperature range of ~ 10 — 10 K that rela- 
tivistic kinetic effects for electrons cannot be ignored. Elec- 
trons diffuse faster than protons (or other nuclei) do as they 
are much lighter in mass, but collectively they are tightly 
trapped by the electric field produced by the difference of 
motions between electrons and positively charged nuclei. 

Matters, including electrons, neutrons, nuclei, are al- 
most completely degenerate inside the highly condensed core 
prior to the onset of a stellar core collapse. However, we con- 
clude from the following calculations [see[Callcn ( 1985)] that 
those materials in the gas shell shortly after a SN explosion 
are far from degenerate. The criterion of non-degeneracy is 
given by (/i is the chemical potential and /3 — (fcaT)"^ with 
fee = 1.3807X 10~^® erg being the Boltzmann constant) 



e^^ « 1, 



or equivalently. 



90 



< 1 



(Dl) 



(D2) 



C alien 



(see 

h = 6.626 X To 



{2nmkBT/h^) ^''^ (here 



19851 where At 

erg s is the Planck constant) is the ther- 



mal de Broglie wavelength, M is the number density of par- 
ticles, m is the particle mass, and po is the intrinsic degree of 
degeneracy; e.g., for particles with a spin quantum number 
s, we have go = 2s + 1. Thus we can derive the degeneracy 
condition of an electron gas (note that M can be estimated 
by A/" — p/(mp^), where rrip is the proton mass and A is 
the average nucleon number per nucleus), viz. 



We find that kinetic motions of electrons are relativistic 
while that of baryons and nuclei being non-relativistic by 
comparing their rest masses and fceT, namely 



m^c = 0.511MeV < k^T < m^c ~ 931MeV < m^ucC^ 



(D4) 



We first estimate the ratio of one nucleus momentum pn 
to one electron momentum pe as 



Pnuc 



(2mnu,fcBT) 



1/2 



(fc|rVc2 



- 7Tl|c^)-'-/^ 



> 



^2m 



1/2 



: 30 



(D5) 



Here k^T roughly gives the energy of thermal motion of one 
particle. The numerical result is based on the assumption 
that the average charge number per nucleus is 2. This es- 
timate implies that we can reasonably ignore the diffusion 
effect of electrons when discussing kinetic effect of diffusion: 
electrons are trapped by the electric field when they go in- 
wards, but they cannot effectively drag the nuclei in since 
Pnuc/Pc 3> 1. Instead, they are dragged by those nuclei; that 
is, we should discuss Pnuc/{Zpc) where Z is the average of 
nuclear charge number, since one nucleus is "bound" with 
Z electrons on average, yet Pnuc/ (Zpc) is still considerably 
greater than 1 for a not- very-large Z. Note that the rebound 
shock has already dissociated a large portion of heavy nuclei 
(e.g. Janka fc Miiller|1996 1. In other words, diffusion effects 
on kinetic aspects are prevailed by the diffusion process of 
nuclei rather than electrons. A similar estimation based on 
the electrostatic force and radiation force is presented by 



e.g. Padmanabhan (20011, yielding similar results 



With a local Cartesian coordinate system erected, with 
the a;— axis pointing radially outwards, at the contact dis- 
continuity interface, probability densities of nuclei velocity 
distribution takes the M-B form of 



p(v) oc exp 



fed) 



2fcRT 



(D6) 



where Vcd is the radial velocity of the contact discontinuity 
interface and m is the average mass of one nucleus. Similar 
to Lou & Zhai (20091, ratio i? of the number of particles 



that diffuses into ro from the region ro < r < ro + 1 {I is the 
mean free path length) during the time period St to the total 
number of particles in that region is approximately given by 



f2Tvm,kBT 
> 2mD -t; 



3/2 



inl3 -3 

10 g cm 



r 



1010 K 



3/2 



(D3) 



where rric is the electron mass. Here we assume that 
the numbers of protons and neutrons in nuclei in the 
gas are roughly the same. This shows that materials in 
the core with as high mass density as p ~ 10^^ g cm""^ 



(e.g. Arnett|[l977 1 are strongly degenerate. Meanwhile, this 
justifies the adoption of the Maxwell-Boltzmann (M-B) 
statistics (i.e. non-degenerate) for electrons at mass den- 
sity p < 10^" g cm~^ ^ 10^'^ g cm~^ and temperature 
T < 10^° K, which are typical values for the innermost re- 
gion of a gas shell being ejected. Similar considerations en- 
able us to apply M-B statistics for baryons and nuclei, which 
are obviously farther from degeneracy since mnuc 2> rric, 
where rrinuc is the average mass per nucleon. 



/ r'>+^ 2 f 3 ^ ^ 

i9 = I / iirr p{r)dr / p(v)d v 

rro+l r 

X / 47rr^p(r)dr / p(v)d^v . 

Jro J |r + t;^<5t|<T'o,(uJ+„2)<(ro/<St)2 

(D7) 

The mean value theorem for integrals enables us to treat p 
and in the integrand as approximate constants, especially 
for those cases with l/ro <^ 1 [this is guaranteed by estimates 
presented by eq. ( |D14[ )], which directly leads to 



1 — exp — 



2fcBT{«t)2 I ro 



ri/2 



rl + l/ro _^ /■[7T(l-*)-"cd] -^ai^T _-2 

da; / e " dt) 



/■l+l/ro r 



(D8) 



[-W(l + «)-"od]\/^ 
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where x = r/ro and v — Vx['m/(2kBT)]^^^ are dimensionless 
integration variables converted from r and v^, respectively. 

One sensible choice of St is 5t — ro/vcd, reflecting the 
time scale during which the void radius almost doubles (t;cd 
is time-dependent in polytropic cases, yet this is still an 
indication of radial flow speed for a period of time). With 
this choice, we have 



nl+ro/i 

X / di [erf(2V + 2x) - erf(Vi)] , 



(D9) 



where we define the dimensionless expansion speed V of the 
contact discontinuity interface 

1/2 



2kBT 



(DIO) 



and erf (a;) is the standard error function of argument x. 

For the scenario of SN remnants, the mean free path I 
of particles near the contact discontinuity interface is 
1 mo 



(Dll) 



nE 

Here E can be either the cross section of particles in 
Coulomb scattering Ecou (scattered by other charged parti- 
cles near the contact discontinuity interface; subscript "cou" 
for "Coulomb") corresponding to Zcou, or Thomson scatter- 
ing Eph (scattered by photons inside the bubble; subscript 
"ph" for "photon") corresponding to Zph. These two kinds of 
scatterings are dominant when evaluating the diffusion pro- 
cesses of matters near the contact discontinuity interface, 
i.e., particles of those matters can be scattered by either 
other charged particles near the shell inner boundary or the 
photons inside the bubble. 

Ecou is estimated in Jackson ( 1999 1 by taking that the 
radius of charge distribution as ~ 1.2 f m x A^''^ 



Ecou ~ 30 f m^Z 



2/3 



hv 



(D12) 



and here e is the unit electric charge, Z is the average 
charge number per nucleon (assuming A/Z ~ 2 and thus 
m ~ 2Zmp), V is to be estimated by the magnitude of 
thermal velocity {v ~ [kBT/{Zmp)Y^'), and h = h/{2Tv) 
is the reduced Planck constant. For T < lO^'' K and Z ^ 1, 
we get Ecou > 10~^^ cm^ (this may be even higher since 
Ecou oc Z^'^^'^) and according to eq. (Dill, Zcou < 10^^ cm. 
For the Thomson case of scattering with photons, we 



have from eq. (51 1 (for the electron-photon scattering) and 



the number density of photons given by Bose-Einstein statis- 



tics (e.g. |Callen||1985[ ): 



in31 -3 

10 cm 



1010 K 



(D13) 



then the mean free path length Iph is also estimated to be 
'ph ~ 10~^ cm as well for the radiation in the bubble whose 
temperature is Trad ~ 10^" K. 

Let us consider a case with these features common for 
SN ejecta when a SN explosion takes place, e.g. ro ~ 10* cm, 
Vcd ~ 10® cm s~^, m ~ 4mp, Ted ~ 10^^° K, and pcd ~ 
These parameters yield following conditions: 

'cou 'd 



10^° g cm"^ 



V ~ 1.6 



tph 
'"0 ''O 



< 10" 



< 1 



(D14) 
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Figure El. A direct comparison between relations of x^d (where 
cicd a-nd 'Jed can be determined) versus xq on the SCC for n = 0.85 
(Panel A) and for n = 0.79 (Panel B). The light solid curves 
marked along with crosses ("x") represent branches of type-1 
void solutions. In Panel A, type-1 void solutions exist for three 
intervals of xq: 0.0094 < xq < 0.47, 1.96 < xq < 2.60, and 2.82 < 
xo < 6.16. In Panel B in contrast, there exists only one branch of 
type-1 void solution for xq in the interval 3.03 < xq < 7.98. We 
find by extensive numerical explorations that there are actually 
no type-2 void solutions when n < 0.87, while the two branches 
of typc-1 void solution on the left and in the middle of Panel A 
also disappear when n < 0.80. 



thus during a time interval of 5t = ro/ucd, 



(1 



erf(2V + 2) - erf(V)] 



1.3% 



(D15) 



This is such a low ratio of paticle diffusion that our void 
model with sharp discontinuous interface is well justified in 
the context of SNe. 

Although this estimation is carried out during a SN ex- 
plosion, we expect it to be hold over a relatively long period 
of time. Since V oc Vcd/T and T oc p/ p for an ideal gas, we 
then have V oc according to eq. ([5|. Also, initial values of 
Ic/ro and Iph/ro are so small that they will not increase to 
be comparable with 1 quickly. 



20 Yu-Qing Lou and Lile Wang 



APPENDIX E: GRADUAL TRANSITION OF 
PHASE DIAGRAM FOR SOLUTIONS 
CROSSING THE SCC SMOOTHLY 

Readers might be interested in the gradual variation of phase 
diagrams indicating the transition from Fig. [5] to Fig. [3] 
e.g., where do some curves go in Fig. |5] as the value of n 
is decreased? In Appendix E here, we briefly show with the 
help of Fig. |E1[ how Fig. [S] becomes Fig. [3] as n decreases. 

As n drops from 0.9, the line marked with circles ("o") 
in Fig. [5j corresponding to the branch of type-2 solutions, 
gradually "shrinks" (i.e. the range of xq for this branch of 
solutions becomes narrower). When n ~ 0.87, this branch 
completely disappears from the phase diagram. Panel A in 
Fig. |E1| shows the phase diagram for n — 0.85 when the 
branch of type-2 solution has already disappeared. 

When n value is reduced further to n ~ 0.805, the left 
two branches of type-1 solutions, indicated by curves marked 
with crosses ("x") in the phase diagram, disappear almost 
simultaneously. That is, only the right-most curve in Panel 
A of Fig. |E1| "survives" . Panel B in Fig. |El| shows the phase 
diagram for n = 0.79 when the left two curves no longer 
exist. Qualitatively, there are no signiflcant changes as n 
continues down to 0.67, except that the range of xq of this 
branch of solutions expands to the extent shown in Fig. [3] 
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